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Application of Maximum Entropy Method to Lattice Field Theory with a 
Topological Term*^ 

Masahiro Imachi a , Yasuhiko Shinno b and Hiroshi Yoneyama c 
a Department of Physics, Yamagata University, Yamagata, Japan 
b Graduate School of Science and Engineering, Saga University, Saga, Japan 
c Department of Physics, Saga University, Saga, Japan 

In Monte Carlo simulation, lattice field theory with a 9 term suffers from the sign problem. This problem can 
be circumvented by Fourier-transforming the topological charge distribution P(Q). Although this strategy works 
well for small lattice volume, effect of errors of P(Q) becomes serious with increasing volume and prevents one 
from studying the phase structure. This is called flattening. As an alternative approach, we apply the maximum 
entropy method (MEM) to the Gaussian P(Q). It is found that the flattening could be much improved by use of 
the MEM. 



1. INTRODUCTION 

It is well known that non-perturbative proper- 
ties of the strong interaction are relevant to the 
dynamics at low energy such as the U(l) prob- 
lem. Although a 9 term is deeply associated with 
the non-perturbative properties, it is indicated 
experimentally that the effect of the 9 term is 
suppressed. This is the strong CP problem. The 
existence of the 9 term also opens the possibility 
of rich phase structures in 9 space. So it is im- 
portant to study the dynamics of QCD with the 
term. 

The 9 term makes Boltzmann weight complex 
in the euclidean path integral formalism. This 
makes it difficult to perform Monte Carlo sim- 
ulation. This problem can be circumvented by 
Fourier-transforming the topological charge dis- 
tribution P(Q). The partition function 2(6) is 
given as 



Q 

where P(Q) is given as 



P{Q) 



S(4>) 



(1) 



(2) 
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The measure [D^q represents that the integral is 
restricted to configurations of the field <j> with the 
topological charge Q, and S denotes an action. 

Although this algorithm works well for small 
lattice volume[|UEl|Sll3, the effect of errors of 
P(Q) becomes serious and disturbs the behav- 
ior of 2(9) as volume increases. In fact, a fic- 
titious signal of a phase transition was observed 
due to the errors of P(Q)[E1 E|- This is called 
flattening. The flattening could be remedied by 
increasing statistics, but it is hopeless because ex- 
ponentially increasing statistics are needed with 
increasing volume. So we do need an alternative 
way El to calculate 2(9) properly. 

In this talk, we use the inverse Fourier trans- 
form and apply the maximum entropy method 
(MEM). As a model, we employ the Gaussian 

P(Q). 

2. MODEL AND FLATTENING 

The Gaussian P(Q) is parametrized as follows; 
P(Q) oc e~^ Q2 . (c, V : parameters) (3) 

This model realizes in the 2-d U(l) gauge model, 
the strong coupling limit of the CP^ -1 model and 
so on. The parameters c and V are regarded as a 
constant depending on the coupling constant and 



2 



a volume in the U(l) gauge model, respectively. 
For our analysis, we use mock data by adding the 
Gaussian noise with the variance of 8 x P(Q) to 
the Gaussian P(Q). 

Fig. ^ shows the free energy density, f(0) = 
— ylogZ(0), obtained by Fourier-transforming 
numerically the mock data for various volumes. 
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Figure 1. Free energy density f(0) for various 
volumes, c = 7.42 and 8 = 1/400. 

As V increases, already at V = 30, f(0) cannot 
be calculated correctly because of errors of P(Q) 
in Fig.[IJ Especially at V = 50, f(0) becomes flat 
for > 2.3 and gives a fictitious signal of a first- 
order phase transition at ~ 2.3. This is nothing 
but the flattening. Since the flattening is char- 
acteristic of the Fourier transform procedure, an 
alternative way should be employed to calculate 
f(0) properly. 



{P(Q)} and information / are given. The infor- 
mation / represents our state of knowledge about 
Z{0). In this case, we impose a criterion that 
Z{0) > 0. 

The probability prob(,Z(0)|P(Q), I) is written 
in terms of \ 2 an d the entropy 5*; 

prob(Z(6)\P(Q),I) a exp{-±x 2 + a S} = e w ^ z \ (4) 

where a is a real-positive parameter. Convention- 
ally the Shannon- Jaynes entropy S is employed [ 



s = 



(W 



Z{9) - m{9) - 2(0) log 



Z{0) 



m{0) 



(5) 



where m(0) is called default model which reflects 
our state of knowledge about Z(0). Namely, 
our task is to explore the image Z{6) such that 
prob(Z(9)\P(Q), I) is maximized by following the 
three steps [ |51 E]; 



1. Maximizing W(Z) for a given a: 
8W{Z) 



8Z{0) 



0. 



(G) 



2. Averaging Z^ a \0) to calculate the 
best image Z(0) which maximizes 
W ob{Z{0)\P{Q)J): 

Z{9)~ I da prob(a|P(Q), I)Z (a) (0), (7) 



3. Error estimation. 



3. MAXIMUM ENTROPY METHOD 

In this talk, the inverse Fourier transform is 
used in the analysis. In this case, the number of 
degrees of freedom of Q, N q , is smaller than that 
in space, Ng, and the analysis by use of the \ 2 ~ 
fit suffers from the ill-posed problem. In order to 
circumvent this problem, we use the MEM, which 
is effective for such the issue. 

The MEM is based on the Bayes' theorem in 
the probability theory. For our case, probabil- 
ity prob(Z(0)|P(Q), I) is considered, which is the 
probability that Z(0) is realized when data of 



4. RESULTS 

In our analysis, the Gaussian P(Q) is used as 
a test. It is a good laboratory to investigate 
how the flattening is improved, because the cor- 
responding partition function Z po i S (0) is calcu- 
lated by use of the Poisson's sum formula ana- 
lytically. We perform the analysis for various pa- 
rameters c, V, 8. Here we fix c = 7.42, V — 50 
and 8 = 1/400. The number of data set is 
30. Three default models are employed: (i) con- 
stant type; m(0) — 1.0, (ii) strong coupling re- 
gion type; m(0) — (|sin|) y , (iii) Gaussian type; 
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m(ff) = exp{— O | 2 10 7^ 2 }7 where the parameter 7 
is varied from 4 to 8. In the analysis, it is non- 
trivial to find a solution due to N„ < Ng and the 
singular value decomposition is employed. In or- 
der to calculate image Z(9) with high precision, 
the Newton method is used with quadruple pre- 
cision. 

Firstly, in order to find a solution which ap- 
proximately agrees with the exact Z po i s (9), analy- 
ses are performed by using the three default mod- 
els for various a by step 1. 
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Figure 2. Z^ a \9) for data at V = 50 with flat- 
tening. As a comparison, the result of Fourier 
transform(o) and Z po is(0)(*) are also displayed. 



Fig. |21 shows the results. Note that all the re- 
sults of the MEM are free from the flattening. 
Results of the constant type at a — 300 and the 
Gaussian one with 7 = 6 at a — 2000 approxi- 
mately agree with the exact Zp 0ls (9). 

In order to investigate whether these solutions 
are favored probabilistically, we calculate Z(9) 
following step 2 and also estimate these errors 
by step 3. We find that Z(9) given by the Gaus- 
sian default model with 7 = 6 is the most favorite 
image in the current analysis. 

The result is shown in Fig. [31 The result is free 
from the flattening and consistent with Z po i s (9). 
Hence we conclude that the MEM works for the 
9 term and reproduces the reasonable image [ITU]. 



Figure 3. Z{9) given by the Gaussian type with 
7 = 6. 

5. SUMMARY 

In this talk, we applied the MEM to mock data 
of the Gaussian P(Q) and showed that the flat- 
tening is much improved by use of the MEM. 

The next task is to apply the MEM to more 
realistic models. 
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